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ABSTRACT 

High-cadence optical observations of an H-a blue- wing bright point near solar AR NOAA 10794 are presented. The data were obtained with the 
Dunn Solar Telescope at the National Solar Observatory/Sacramento Peak using a newly developed camera system, the RAPID DUAL IMAGER. 
Wavelet analysis is undertaken to search for intensity-related oscillatory signatures, and periodicities ranging from 15 to 370 s are found with 
significance levels exceeding 95%. During two separate microflaring events, oscillation sites surrounding the bright point are observed to twist. 
We relate the twisting of the oscillation sites to the twisting of physical flux tubes, thus giving rise to reconnection phenomena. We derive an 
average twist velocity of 8.1 km/s and detect a peak in the emitted flux between twist angles of 180° and 230°. 
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1. Introduction 

Magnetic reconnection remains one of the most promising 
mechanisms to convert magnetic energy into heating of the 
local plasma (see Klimchuk 2006 ). Indee d, recent theoretical 
work by De Moortel & Galsgaard (i2006ah suggests that photo- 
spheric footpoint motions lead to the build-up, and subsequent 
release, of magnetic energy. Two different types of footpoint 
motion have been proposed as potential mechanisms for re- 
connection, large-scale rotational motion and small-scale spin- 
ning. The proposed mechanism relies on two initially aligned, 
thin flux tubes, which are forced to interact due to foot- 
point motions. Twisting of the magnetic field in such a sys- 
tem will lead to the linear coupling of various MHD modes. 
However, the degree to which magnetic twisting modifies the 
dispersion relationship is not well understood (Nakariakov & 
Verwichte 2005). 



Recent work by Leenaarts et al. (I2005h has discussed 
the formation of the H-a wing in detail using 3D magneto- 
convection simulations, and they show that H-a blue-wing 
bright points correspond to kiloGauss magnetic field concen- 
trations in the photosphere. Co-aligned hard X-ray, H-a and 
magnetic field observations show that the X-ray emission of 
microflares occurs in small magnetic loops, while H-a emis- 
sion originates from the footpoints of these loops located in 
the lower solar atmosphere (Liu et al. 2004). Due to the coinci- 
dence of these bright points with high magnetic field concentra- 
tions, magnetohydrodynamic (MHD) waves (Kalkofen 1997; 
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Fig. 1. The H-a blue wing field-of-view obtained at the DST in solar 
heliocentric coordinates. The colour scale indicates normalized flux 
and the solid box outlines the 40 x 50 pixel region removed for sub- 
sequent bright point analysis. The two dashed boxes outline additional 
bright points which were analysed for control purposes. 



Hasan & Kalkofen 1999) appear to be the modes of wave prop- 
agation likely to be present (McAteer et al. 2002). Theory sug- 
gests that forms of MHD wave with longitudinal components 
to their wave- vector can be induced via reconnection events 
(Hollweg 1981). 

If reconnection processes are to be both temporally and 
spatially resolved, it is imperative to acquire data at the high- 
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Fig. 2. The AR bright point lightcurve where the intensity is com- 
puted by averaging pixel values which supersede a lower-intensity 
threshold of the background modal value plus 10a. The vertical 
dashed line at 630 s corresponds to the time when the oscillatory 
power has undergone a 180° twist. The second vertical dashed line 
at 685 s corresponds to the time when the oscillatory power has un- 
dergone a 230° twist. Notice how the peak flux falls between the two 
vertical dashed lines in accordance with the theory in § |4] 




10 20 30 10 20 30 

Pixel Number Pixel Number 



est possible cadence using a sensitive camera system. Here we 
report intensity oscillations originating from an active region 
( AR) bright point, and the corresponding spatial rotation during 
microflare activity detected within the lower solar atmosphere. 

2. Observations 

The data presented here are part of an observing sequence ob- 
tained on 2005 August 10, with the Richard B. Dunn Solar 
Telescope (DST) at Sacramento Peak. Use of a highly sensitive, 
dual camera system enabled synchronized, simultaneous im- 
ages at a rate of 20 frames per second in each camera to be ob- 
tained. The RAPID DUAL IMAGER (RDI) camera system, used 
to acquire the observations presented here, was developed by 
Queen's University Belfast and remains stationed at the DST as 
a common user instrument. The RDI system comprises of two 
Basler A301b cameras - one 'master' and one 'slave' channel 
- connected to a custom built PC, and controlled by software 
designed and developed by 4C's of Somerset, UK. The cam- 
eras have a 502 X 494 pixel^ CCD, with a square pixel size of 
9.9 /im and can operate at a maximum speed of 80 frames per 
second. Currently, the PC has 200 GB of available disk-space 
allowing for an uninterrupted data run of over 7 hours. 

Our optical setup allowed us to image a 50.4'' x 49.2'' re- 
gion surrounding AR NOAA 10794 complete with solar ro- 
tation tracking. The active region under investigation was lo- 
cated at heliocentric co-ordinates (770",-254"), or S12W56 
in the solar NS-EW co-ordinate system. A Zeiss universal bire- 
fringent filter (UBF; Beckers et al. 1975) was used for H-a 
blue- wing (YL-a core - 1.3 A) imaging with one of the RDI 
CCD detectors. The 6561.51 A H-a blue- wing central wave- 
length selected with the UBF had an associated filter bandpass 
of 0.21 A (Beckers et al. 1975). In addition, a G-band filter was 
employed with the second RDI camera to enable synchronized 
imaging in the two wavelengths. Due to the smaller granulation 
contrast of H-a blue- wing images when compared to G-band 
data, only results obtained in the H-a blue wing are presented 
here. Leenaarts et al. (2005) have demonstrated that, as a con- 
sequence of this reduction in granulation contrast for H-a blue- 



Fig. 3. The AR bright point before (upper left) and during (upper 
right) the first microflaring event. The lower left panel displays the 
bright point prior to the second microflaring event, while the lower 
right panel shows the bright point undergoing the second microflaring 
event. Overplotted in each diagram (green contours) are locations of 
strong (> 80% of the maximum of the entire time series) oscillatory 
power at a periodicity of 20 s. The microflaring bright point images 
reveal strong oscillatory power surrounding their centre as described 
in §131 



wing images, better feature tracking is obtained. During the ob- 
servational sequence used in the analysis, low-order adaptive 
optics was implemented. 

The data selected for the present analysis consist of 31760 
H-a blue- wing images taken with a 0.05 s cadence over a total 
time period of 26.5 min. The acquisition time for this observ- 
ing sequence was early in the morning and seeing levels were 
good with minimal variation throughout the time series. These 
images have a spatial sampling of 0.1" per pixel, to match the 
telescope's diffraction limited resolution in the H-a blue wing 
to that of the CCD. Further detail s rega rding the RDI camera 
system can be found in Jess et al. torn . 



3. Data Analysis 

Small-scale turbulent seeing in the Earth's atmosphere means 
that even high-order adaptive optics cannot compensate for all 
rapid air motions. Here we implement the speckle reconstruc- 
tion masking method of Weigelt & Wirnitzer (1983), adapted 
for solar im aging b y von der Liihe (1993) and further improved 
by de Boer ( 1995k . Observing at a high cadence, the short ex- 
posure times essentially freeze out atmospheric distortions and 
maintain signals at high spatial frequencies, albeit with statis- 
tically disturbed phases (Siitterlin et al. 2001). It is possible 
to recover the true amplitudes and phases in Fourier space by 
taking a large number of such short exposure images, called a 
"Speckle Burst", and utilizing an elaborate statistical model. 

Eighty raw RDI data frames were used for each speckle re- 
construction, producing a new effective cadence of 4 s. This 
provides a Nyquist frequency of 125 mHz and is therefore 
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Fig. 4. The normalized flux of the bright point during the second microflare event is plotted, as a function of time, in the top panel. The bottom 
panel shows the spatial representation of oscillatory power for 20 s periodicities. The letters placed above each oscillatory power map indicate 
the time, through comparison with the top panel, at which the spatial wavelet-power image is generated. It can be seen that the twisting of the 
oscillation sites is closely linked to the brightening of the bright point. This figure is available as a movie for the online version. 



suitable for the search of oscillations with periods longer than 
8 s. Typical Fried parameters obtained prior to speckle recon- 
struction were ro ~ 10 cm, indicating good post- speckle im- 
age quality. Mikurda & von der Liihe (20061) have shown that 
although the residual wavefront errors differ between images 
obtained using low-order adaptive optics and uncompensated 
data, its effects on the speckle transfer function are, to a first- 
order approximation, similar to that of an uncompensated wave 
obtained during better seeing. Thus, since speckle reconstruc- 
tion takes all residual aberrations into account, implementa- 
tion of the speckle masking method outlined above is possi- 
ble, on data obtained using low-order adaptive optics, without 
additional processing. However, after implementing temporal 
Fourier analysis on 2400 successive dark images, a moving- 
pattern noise was revealed, showing up as an intense Fourier 
power with a variable periodicity between 3 and 9 s. Therefore, 
to strengthen the reliability of our detections, only oscillations 



with periodicities greater than 15 s were incorporated into the 
analysis. Before commencing Fourier and wavelet techniques, 
a 40 X 50 pixel (4" xb") region surrounding the AR bright point 
under consideration (Fig.[T]) was isolated for further study. All 
other spatial information was discarded to allow the analysis to 
focus on the behaviour of the bright point. 

To compensate for camera jitter and large-scale air-pocket 
motions, all data was subjected to a Fourier co-aligning rou- 
tine. This routine utilizes cross-correlation techniques as well 
as squared mean absolute deviations to provide sub-pixel co- 
alignment accuracy. However, it must be noted that sub-pixel 
image shifting was not implemented due to the substantial 
interpolation errors which may accompany use of this tech- 
nique. Since the observing sequence was obtained in the early 
hours of the morning, when image warping is particularly 
strong, all data were de-stretched relative to simultaneous, 
high-contrast G-band images. We use a 40 x 40 grid, equat- 
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Fig. 5. A trace of the oscillation sites after de-rotation to solar disk 
centre. The path taken is not circular, providing a maximum traversed 
distance of 1230 km. 



ing to a 1.25^' separation between spatial samples, to evaluate 
local offsets between successive G-band images. Due to both 
cameras sharing the same pre-filter optical path, all determined 
local offsets are applied to simultaneous narrowband images 
to compensate for spatial distortions caused by atmospheric 
turbulence and/or air bubbles crossing the entrance aperture 
of the telescope. However, even though wave front deforma- 
tions in the pupil plane are first order wavelength indepen- 
dent, the image degradations are not, as can be seen when 
comparing derived Fried parameters at different wavelengths. 
Nevertheless, the de-stretching grid implemented in this pro- 
cess allows moderate- to large-scale seeing-induced distortions 
to be removed, since they are less wavelength dependent. With 
this in mind, very small-scale atmospheric distortions may not 
be accurately compensated for. 

After successful co-alignment and de- stretching, time se- 
ries were created for each pixel before being passed into Fast 
Fourier Transform (FFT) and wavelet analysis routines. The 
wavelet chosen for this study is known as a Morlet wavelet, 
and is the modulation of a sinusoid by a Gaussian enve- 
lope (Torrence & Compo 1998). A number of strict crite- 
ria implemented on the data allowed us to insure that oscil- 
latory signatures correspond to real periodicities. These cri- 
teria have been described in detail in previous papers (see 
Banerjee et al. 2001, McAteer et al. 2004, Ireland et al. 1999, 
Mathioudakis et al. 2003). 

Four-dimensional maps containing spatial information as 
well as wavelet power and oscillatory period were saved as 
outputs of wavelet analysis for detected oscillations which un- 
der the above criteria had a greater than 95% significance 
level. Furthermore, a time series was created for the AR bright 
point to display its temporal variability throughout the time se- 
quence. The bright point was identified by establishing which 
pixels in the 40 x 50 pixel region displayed intensity levels ex- 
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Fig. 6. The rotational angle of the footpoints as a function of time. 
The triangular and star symbols indicate the nature of the twisting for 
each of the oscillation sites, while the bold solid line shows the average 
rotational angle. As described in § |4l there is clearly an acceleration 
of the twist velocity between an angle of 20 and 160 degrees. Above 
160 degrees, the twist velocity decelerates continually until the motion 
ceases. 



ceeding the background median value plus lOa. In each frame, 
approximately 30 pixels exceeded this threshold, indicating a 
bright point surface area greater than 150000 km^. The pix- 
els which superceded this intensity thresholding were averaged 
and placed at the corresponding time element in the time se- 
ries array. Hence, variability in oscillatory phenomena can be 
compared readily with the bright point flux. 

4. Results and Discussion 

An examination of the bright point time series over the com- 
plete 26.5 min duration of the observations reveals two events 
which show intensity increases (Fig.O. The first corresponds 
to an approximate 80% intensity increase over the neighbour- 
ing quiescent flux, with an event duration exceeding 120 s, 
while the second is a much larger event and relates to an 
approximate 125% increase in intensity lasting in excess of 
350 s. Leenaarts et al. (2006) have shown that the brightness 
of H-a blue-wing features, just as for the G-band, will de- 
pend on the orientation of flux tubes, with higher intensities 
seen when the observer is looking straight down the tube axis. 
Thus, the swaying of flux tube bundles might explain tempo- 
rary brightness enhancements when the buffeting of granules 
pushes the flux tubes, which are normally pointing 50 degrees 
away from the observer, towards them. However, due to the 
sharp rise in intensity and the corresponding decay pattern, 
we deem these intensity fluctuations to be characteristic of mi- 
croflare activity with similar energetics and durations observed 
by Porter et al. (Il987k . 

A succession of wavelet power diagrams around the time at 
which the microflaring activity occurs reveals a periodic signal. 
Prior to the commencement of both microflare events, there is 
an increase in oscillatory power superimposed, spatially, over 
the bright point (Fig. [3]). The increase in oscillatory power pro- 
vides two distinct regions of power separated by the central 
portion of the bright point. This increase in oscillatory power 
occurs 120 s prior to the second, larger microflare event which 
corresponds to 6 complete oscillation cycles for the 20 s pe- 
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Fig. 7. Phase difference between the two oscillating sites, regardless 
of the magnitude of oscillatory power, prior, during and after the twist- 
ing is seen. Note the lack of anti-phase oscillations and minimal phase 
variation throughout the time series, particularly during the time inter- 
val when the second microflare event occurs (500-850 s). 



riodicity plotted in Figure [3] During the evolution of both mi- 
croflares, a spatial twist of the oscillatory sites occurs. The os- 
cillation sites are not symmetric, yet appear to pivot around the 
bright point centre. Focusing on the second microflare event, 
Figure |4] reveals the nature of the twisting along with the asso- 
ciated time series. 

Due to the position of the active region near the West 
limb, line of sight effects must be considered. The observed 
circular path taken by the rotating oscillation sites will actu- 
ally correspond to an elliptical path after consideration of the 
cosO = 0.62 term. The path taken by the oscillating sites, after 
de-rotation to disk center, is shown in Figure [5] The perime- 
ter distance around this path is 17 pixels, corresponding to a 
traversed distance of 1230 km, while the duration of twisting 
is 148 s. The rotational movement of the oscillatory twist fails 
to complete a full 360° rotation, but we can use this distance 
to derive an upper limit to the average rotational velocity of 
the twisting oscillation sites, i.e. the twist velocity. Simply di- 
viding the maximum traversed distance (1230 km) by the time 
taken (148 s), we establish a maximum average twist veloc- 
ity of 8 km/s, very close to the photopsheric sound speed 
(Fossum & Carlsson 2005). To demonstrate the temporal char- 
acteristics of the twist velocity. Figure [6] plots the rotational 
twist angle as a function of time. From this figure it is clear 
that the rotational motion is accelerated between twist angles 
of 20° and 160°, before being decelerated at rotational angles 
exceeding 160°. The twist angle is generated by tracking the 
wavelet power maxima as the oscillation sites pivot around the 
centre of the bright point. By establishing the movement of suc- 
cessive oscillatory power maximas, relative to the central pivot 
point, we are able to compute the twist angle directly. The high 
cadence of the observations, coupled with the high oscillatory 
power, provides a direct, and reliable, measure of the twist an- 
gle. 
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Fig. 8. The temporal evolution of the bright point (image scale in pix- 
els). The colour scale indicates the number of times a pixel has disen- 
tangled itself from the confines of the bright point perimeter. The thin 
black contour outlines the average position and shape of the bright 
point during the entire observing sequence, while the thick black con- 
tour shows the time-averaged position of the oscillation sites. 



We have carried out a number of rigorous tests to assess the 
reliability of the observed twisting. The first test investigates 
the phase relationship of the two oscillating sites by compar- 
ing the positions of intensity peaks and troughs through use of 
a cross-power spectrum. Had the results indicated anti-phase 
behaviour, then the origin of the opposing oscillations would 
have been dubious. However, from Figure [71 it is clear that 
minimal anti-phase exists between the two oscillating sites, in- 
dicating real oscillatory phenomena. This finding is intuitive, 
since if the oscillation sites mark the footpoints of the same 
flux tube, any impulsive energy deposited, through reconnec- 
tion, would produce oscillations which are coupled in phase 
(Nakariakov et al. 1999, Wang et al. 2003). 

The second test creates a mask which maps the evolution of 
the bright point. Since intense pixels may move in and out of 
what we deem to be the perimeter of the bright point, this may 
produce large oscillatory power around its confines. We cre- 
ate a bright point pixel-mask for each frame of the time series, 
and by examining the image differences to previous masks, we 
are able to study the long-term evolution of the bright point. 
By summing all positive movements detected from successive 
masks, i.e. outflows, we can generate a picture of when, where, 
and by how much the bright point shape changes (Fig. [8]). After 
examination of Figure [H we can see some long-term evolution 
of the bright point, particularly around the North-Easterly edge 
where 7 outflows are detected over the entire 26.5 min dura- 
tion of the time series. Overplotting the times related to these 
bright point outflows, on top of a wavelet power transform for 
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Fig. 9. Wavelet transform for the region of the bright point which undergoes the most temporal evolution during the 26.5 min time series. The 
original pixel lightcurve, normalized to zero, is shown in a), with the wavelet power transform along with locations where detected power is at, 
or above, the 99% confidence level are contained within the contours in b). Plot c) shows the summation of the wavelet power transform over 
time (full line) and the Fast Fourier power spectrum (crosses) over time, plotted as a function of period. The global wavelet (dotted line) and 
Fourier (dashed dotted line) 95% significance levels are also plotted. The cone of influence (COI), cross-hatched area in the plot, defines an 
area in the wavelet diagram where edge effects become important and as such any frequencies outside the COI are disregarded. Vertical lines 
indicate the times corresponding to outflows from the bright point. Notice how the outflows do not coincide, temporally, with the creation of 
high 20 s oscillatory power, which correspond to the micro-flare events described in § |4] 



the corresponding flux of the North-Easterly region, reveals no 
correlation with the locations of intense 20 s oscillatory power 
(Fig.O. Furthermore, on closer inspection of Figure [SI we see 
that the time-averaged position of the oscillation sites occupy 
a larger area than the bright point, indicating no dependence of 
the oscillation sites on outflows and/or jitter of the bright point. 

Finally, our third test involves the examination of the bright 
point flux. Even though we see the oscillation sites pivot around 
the bright point centre, some of these oscillation sites overlap 
with its confines. Thus, we should be able to detect intensity 
oscillations, albeit with much reduced power due to averag- 
ing over bright point pixels which are not oscillating, in the 
normalized flux output of the bright point. Running the bright 
point time series through wavelet analysis techniques, using the 
same strict criteria outlined in § [S] we see that the 20 s oscil- 
lation is detectable (Fig.[TO|. This enforces our belief that the 
detected oscillations are real, and we are not seeing an artifact 
of bright point evolution and/or jitter. 

Leenaarts et al. (2006 ) have shown that the H-a wing is a 
great tool to locate and follow long-lived isolated magnetic el- 



ements, and since bright points in the blue wing of the H-a line 
correspond to kiloGauss magnetic field concentrations, mag- 
netic features, such as flux tubes, may be present. Introducing 
the theoretical work by De Moortel & Galsgaard (2006a), we 
can interpret the oscillatory sites described above as a bundle 
of initially aligned, thin flux tubes. Even though we cannot di- 
rectly see the flux tubes, material flowing along such paths may 
induce magneto-acoustic oscillations, which we can detect with 
the wavelet analysis applied to the data. The De Moortel & 
Galsgaard ( 2006b) model relates the evolution of the flux as- 
sociated with magnetic reconnection to the rotational angle of 
initially aligned flux tubes. 

Interpreting our results as having a non-zero background 
potential field (consistent with same-polarity flux domains), the 
peak reconnection rate should occur between a twist angle of 
180° and 230°. Figure [2l shows the bright point time series dur- 
ing the second microflare event, along with vertical lines indi- 
cating the times where the twist angle equals 180° and 230°. 
From this plot it is clear there is a peak in the emissive flux 
between these two angles. On closer examination of Figure [6l 
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we can see that during peak acceleration, the twist velocity 
is greater than the sound speed. Between a time of 560 and 
600 s, the oscillation sites are pivoted by approximately 160°, 
corresponding to a traversed distance of 540 km, and provid- 
ing a velocity approaching 14 km/s. Kondrashov et al. d 1999b 
have studied three-dimensional magnetohydrodynamic simu- 
lations of magnetic flux tubes, and have concluded that foot- 
points of interacting flux tubes, which are allowed to move to 
simulate conditions in the solar photosphere, promote the cre- 
ation of supersonic motions during reconnection. We indeed 
see, from Figure O the development of supersonic motion after 
commencement of a micr oflare, which is consistent with the 
work by Kondrashov et al. (|l999|) . 

Examining Figs. [6l and [21 more closely, it appears that the 
impulsive phase of the bright point microflare commences with 
a minimal twist angle. This is not in agreement with the model 
of De Moortel & Galsgaard (l2006bl) . however, explanations can 
be proposed to help understand why th is is the case. The model 
of De Moortel & Galsgaard (l2006bh relies on two initially 
aligned, line-tied plates between which the magnetic loops are 
not curved. This constrains the model with arbitrary parameters 
which may not accurately represent the observations presented 
here. Additionally, the distance of the plates between one an- 
other fully determines the degree of magnetic twist per unit 
length, and coupled with the loop curvature, predicts the start- 
ing time of reconnection. Furthermore, the arbitrary treatment 
of key input parameters such as diffusivity and pressure add 
to potential discrepancies between the model and the obser- 
vations presented here. Thus, the model may not be advanced 
enough to accurately predict the minimal twist angle required 
for reconnection. Contrarily, twisting of the flux tubes may 
have commenced prior to our detection, with this discrepancy 
caused by our RDI camera system having insufficient sensitivity 
to detect such weak signals. Current generation optical detec- 
tors should be able to address this problem with their improved 
dynamic ranges and sensitivities. This would allow weak rota- 
tional signals to be studied and more accurately compared with 
theory. Referring back to the first microflare event, a similar 
progression of oscillatory site evolution occurs. However, this 
event is less energetic, and as such our ability to trace exact 
movements is greatly reduced. 

Two other bright points in the same field of view were also 
analysed using the same criteria described in § [S] These bright 
points are of similar size and luminosity to the bright point in 
question (Fig.[T]). We find no evidence for flaring in these two 
locations and there is no twisting of oscillatory power. 

5. Concluding Remarks 

We have presented direct evidence of high-frequency magneto- 
acoustic oscillations occurring in the immediate vicinity of an 
H-a blue- wing bright point. Periodicities as short as 15 s are 
found with significance levels greater than 95%. We have in- 
terpreted the rotational movement of oscillatory sites during 
microflare activity as the physical signature of flux tube twist- 
ing. We derive the rotational angle corresponding to maximum 
flux and the average twist velocity. This observational result is 
in qualitative agreement with the numerical model introduced 
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Fig. 10. Wavelet transform of the normalized flux of the bright point. 
All scales and axes are identical to Figure [9l except for the inclusion 
of a probability plot in d). The probability of detecting non-periodic 
power is calculated for the peak power at each timestep by comparing 
the value of power found in the input lightcurve with the number of 
times that the power transform of a randomized series produces a peak 
of equal or greater power. A percentage confidence is attributed to the 
peak power at every time step in the wavelet transform such that a low 
value of probability implies that there is no periodic signal in the data, 
while a high value suggests that the detected periodicity is real. We 
see that during the time of the second microflare event, strong power 
with high probability values, is associated with the detection of a 20 s 
periodicity. 



recently by De Moortel & Galsgaard (2006a,b). A direct quan- 
titative comparison with the model is not possible due to uncer- 
tainties associated with key input parameters such as resistivity, 
diffusivity, pressure, the distance between opposing line-tied 
plates and the corresponding degree of magnetic twist per unit 
length. 
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